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Abstract. We define the topological entropy per unit volume in parabolic PDE's such as the complex Ginzburg- 
Landau equation, and show that it exists, and is bounded by the upper Hausdorff dimension times the maximal expansion 
rate. We then give a constructive implementation of a bound on the inertial range of such equations. Using this bound, 
we are able to propose a finite sampling algorithm which allows (in principle) to measure this entropy from experimental 
■ data. 



1. Introduction 



In this paper, we shall deal with a general reaction-diffusion equation, and we have in mind an 
iV-component system in R d which is of the form 

N 

d tUi (x, t) = J2d ij Au j + F i (u v ..., u N ) , i = l,...,N, (1.1) 

where all quantities are real. For example, the complex Ginzburg-Landau equation (CGL) is 

t)+v(x,t) - (l + if3)v(x,t)\v(x,t)\ 2 , (1.2) 

which clearly can be brought to the form of Eq.(l . 1) by writing equations for the 2 components 
u l = Rev, u 2 = Imv. We shall write Eq.(l.l) short as 



d t u = 



DAu + F(u) 
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We state now our assumptions on D and F. We let D sym = (D + D T ) jl denote the symmetric 
part of D, and we assume that the matrix D sym has spectrum {fj}^^...^ in the open right 
half-line. We next define 

v * = . m in M , D * = \\ D \\ , ( L3 ) 

% — 1, ...,iv 

where ||D|| is the norm of the matrix D as a linear map from (equipped with the I 2 norm) 
to itself. Note that is the minimal dissipation rate in Eq.(l.l). Our assumptions on F are 
somewhat vague, but they are intended to cover a large variety of possible applications. We first 
assume that F is "globally stabilizing" in the sense that there is a constant such that for any 
initial condition u(x,0) which is bounded in L°° there is a t < oo such that \\u(-, i')ll<x> < 
for all t' > t. (The factor i is convenient for later use.) In this sense, is the radius of a 
globally invariant set (usually this will be an attracting set). 1 In the case of the real Ginzburg- 
Landau equation (Eq.(1.2) with a = (3 = 0), one has Q*/2 = 1 and for the CGL it has been 
shown in [C,GV] that < oo when d = 1 or d = 2, and also in dimension d = 3 for some 
nontrivial parameter range of a and (3. For many other equations one can derive similar bounds 
using the localization techniques of [CE1]. Since this is not the central issue of our paper, we 
shall just assume that 

KIL<Q*/2, (1.4) 

for all t. Here, and in the sequel u t (x) = u{x,t). Our last general assumption is a bound on the 
maximal local expansion rate. Consider two solutions u and v with \\u t and \\v t bounded 
by for all t > 0. (This is no loss of generality if we consider later functions in the "global 
attractor" A) We define w = u — v. Then, we assume that the non-linearity F is such that w 
satisfies an equation of the form 

d t w(x, t) = DAw(x, t) + M(x, t)w(x, t) , (1.5) 

where the matrix M has a norm (as a map from to itself) bounded by 

||M(M)|| < M., (1.6) 

for all x and t. If F is a polynomial, such a bound will follow automatically from the bound of 
Eq.(1.4). For example, for the CGL, written in complex notation, we have 

d t w = (1 + ia)Aw + (1 + (1 + i(3)(uv + vv))w + ((1 + i/3)u 2 )w , 

so that in this case M is a 2 x 2 matrix whose norm is bounded by 

||M|| < 1 + 3(1 + \P\)Q 2 J4. 



We shall try to stick to the following notation: Quantities with a * as an index depend on the parameters of the 
Eq.(l.l), and the constants C , C v . . . do not depend on them. 
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Again, many other examples can be handled in this manner and are left to the imagination of 
the reader. 

Our study of the topological entropy is based on a detailed analysis of the Eq.(1.5), and in 
particular on the control of information which is localized in space, in the conjugate momentum, 
and in time. The localization in space and time has been developed earlier [CE2] and used to 
prove the existence of the e-entropy per unit volume of Kolmogorov and Tikhomirov [KT] for 
systems such as CGL. Here, we use these estimates to show the existence of the topological 
entropy per unit volume. We then improve the bounds to localize at high frequencies, where 
the flow defined by Eq.(1.5) will be seen to be essentially a contraction. Using this information, 
we shall then show that the topological entropy per unit volume can be measured in terms of 
a discrete sampling of the signal u in space and time. The amount of data needed for such an 
enterprise is, however, quite prohibitive [ER2] if any reasonable precision is to be attained. But 
this is probably unavoidable. For a similar study in finite volume, see [CJT]. 

We end this introduction by explaining in more detail how the various physical scales 
interact as we bound w, since this should be useful to prove further results for dissipative 
systems in unbounded domains. We wish to argue in "dimensionally correct units" so that w 
has the dimension of the observed fields (e.g., a temperature), [M] ~ [t]~ l and [D] ~ [£ 2 t~ l ], 
where t is time and £ is length. In the long wavelength limit, diffusion is inactive, and the growth 
of w is dominated by M. Given the a priori bound on M, we shall fix the unit of time to 

(L7) 

Recall that is the norm of D and the (dimensionless) quotient D^/v^ compares essentially 
the strongest to the weakest dissipation rates. The time 7% is the time in which errors can grow 
at most by a bounded factor, which depends on D^/u^, and we shall see that l/r + is also the 
natural sampling rate for the determination of the entropy. The dissipative range of the equation 
(1.5) is given by those /c-values for which dissipation is guaranteed to exceed the growth, i.e., 
for i/ + k 2 > . Hence we set the cut-off for the fc-values to 

K = (-±-Y /2 -F(DjvJ, (1.8) 

' * * 

where the correcting factor F will be defined in Eq.(6.12). The natural unit of length is almost 
the inverse of k^ : 

6. - (r^Y /2 - (1-9) 

(The correcting factor in Eq.(1.8) is used in the bounds, but for a more intuitive understanding 
the reader should assume = v+.) In terms of these units, we can now explain our "sampling 
bound" of Theorem 6.7: Consider an e > 0, which will be the precision we want to achieve (up 
to a factor). Assume that two solutions of Eq.(l.l), it, and v satisfy the bound 



u(m5 m ,t - nrj - v(m5 m ,t - nrj| < e, 



(1.10) 
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for all n with < n < C \og(l/e) and for all \m\ < (L + C/e)/5^ where C is some universal 
constant. There is another universal constant C' such that Eq.(l.lO) implies 

\u(x,t) - v(x,t)\ < C'e, (1.11) 

for all | x | < L. Thus, a discrete sampling at spacings 5^ and suffices to bound the difference 
of the two functions everywhere in {\x\ < L}. In other words, sampling in the inertial range 
for a time of order log(l/e) suffices to bound the dissipative part of u — v as a function of 
the sampling precision. 




r,B* 5 log(l/e) 



L + Bye 

Fig. 1: A symbolic representation of the space-time points which need to be sampled with differences less than e to 
get a bound on the difference of two functions on the interval labeled L to a precision B 25 s. See Theorem 6.7 for a 
definition of the constants B^, B^ 4 and B^ 5 . 

This result is in line with our earlier bound [CE2] where we showed that, expressed in 
the units of the present paper, the number of balls of radius e in L°° needed to cover the global 
attractor (when restricted to \x\ < L) has a bound of the order of 

exp(C(L/5J d - 1 log(l/e) 2 ) .exp(C(L/5J d log(l/ £ )) . 

One of the main ingredients of the proof of Eq.(l . 1 1) is the "dissipative bound" of Theorem 
5.1 below. In [CE2], we showed the inequality [CE2, Eq.(6.6)]: 



sup 

|x|<A-e- 



\w{x,t)\ < Ce , 



(1.12) 
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when < t < r + , and sup| x | <A \w(x, 0)| < e. We can improve this bound (slightly) as follows 
and write it in natural units (see Eq.(5.8) for a precise formulation): 

sup \w(x,tJ\ < C(sup \w(xM + (K/m^(-MU( — ) {i+d)/2 - 

\x\<\-£ \x\<\ u * 

Our improved bound in this paper exploits the dissipation as follows: Let P > denote a localization 
of the Fourier transform of w to wave numbers k satisfying \k\ > k^ (see below for a more 
precise definition). Then, 

sup \P >W {x^)\ < C(^y~ d /\ sup k(a;,0)| + (^)|K.,0)|| oo )(^)( 1 + d )/ 2 . 

\ x \<x-e * * |x|<a * 

Thus, if k^ is sufficiently large, the high frequency part of w decays. 

Our paper is organized as follows: We first show that the topological entropy per unit 
volume is defined and is bounded by the expansion rate. We then show the dissipative bound 
mentioned above, and finally show how it can be used to prove the sampling bound. From this, 
an estimate of the topological entropy per unit volume can be derived. 

Acknowledgments. This work was partially supported by the Fonds National Suisse. Our 
collaboration was also made possible through the pleasant atmosphere at the IHES, Bures-sur- 
Yvette. 



2. Existence of the Topological Entropy per Unit Volume 

We start by proving existence of the topological entropy per unit volume. This is somewhat 
similar to the standard proof of existence of topological entropy (see [AKM]). We define the 
global attractor A by 

A(B,r) = f)$ nT (B). (2.1) 

n>0 

Here, B is a ball in L°° and $ t is the semi-flow defined by the evolution Eq.(l.l). It can be 
shown (see [MS]) that A(B, r) is invariant and that it does not depend on the initial ball B (if it 
is large enough) nor on the (large enough) time r > r Q (B). Thus, we define A = A [B, r (S)) . 

For any e > and any bounded set Q in R d , whose boundary has Lebesgue measure 0, we 
define Wq as the set of all finite coverings of A by open sets in L°°(Q) of diameter at most e. 
Note that by the compactness of A\q, which follows from the uniform bounds on the gradient 
by Ascoli [MS], such finite coverings exist. 

Let r > be a fixed time step, and let T = m with n E Z. For U e Wq, we say that two 
trajectories A l and A 2 in A are W-different before time T if there is at least one j, < j < n 
for which the points $ JT (Aj) and $>-(A 2 ) do not belong to the same element of U. Let 



N TtT (U) 
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be the largest number of trajectories which are pairwise W-different before time T (and con- 
sidered with time-step r.) Note that this number is finite since it is at most (Card U) T ^ T . 
Let 

N QTre= inf N Tt {U). 



Lemma 2.1. Let Q l and Q 2 be two bounded domains with boundary of zero Lebesgue 
measure and Q l D Q 2 of zero Lebesgue measure. The functions Nq T satisfy the following 
bounds: 

i) Nq Tt e is non-increasing in e. 

H) N Q,T l+ T 2 ,T,e < N Q,T u T,e N Q,T 2 ,r,e- 
m ) N Q l UQ 2 ,T,r,e ^ N Q 1: T,t,s N Q 2 ,T,t, £ - 

Proof, i) follows since Wq is non-decreasing in e. For a given U £ Wq, we have easily from 
the definition 

n Ti+ t 2 ,A u ) ^ N Tl A u ) n t 2 A u ) ■ 

Indeed, if we consider a maximal collection of trajectories A l , A 2 , ■ ■ ■ , A Nt +t we can 
collect with A l all the remaining A - (j > 1) whose step r orbit U coincides with that of A l up 
to Tj . We know that since these orbits should differ between Tj and T 2 their number is at most 
N T ^ T {U). We continue with the remaining orbits and this leads to at most iV Ti T (U) groups, ii) 
now follows immediately. 

In order to prove iii), we consider^ £ Wq and U 2 £ Wq 2 - Since we are using the L°° 
norm, an argument similar to the one above leads to 

N TtT {U,xU 2 ) < N TtT {U x ) N TjT (U 2 ) . 
We also have easily (again because we are using the L°° norm) 

W Q] x Wq 2 C W QiU q 2 . 

The claim iii) now follows easily. 

Definition. Henceforth, we shall work with domains which are cubes, and we denote Q L the 
cube of side L centered at the origin. 

Theorem 2.2. The following limit exists 

h = lim lim — T lim 1= log N n rp . (2.2) 

Moreover, h does not depend on t. It is called the topological entropy per unit volume of the 
system. 

Proof. From ii) of the Lemma 2.1, and the usual subadditivity argument, we conclude that the 
following limit exists 

hn\-r = lim -logiV nTTf . 
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Moreover, from i) it is non-increasing in e and from iii) it satisfies 

h {1) < h (1) +h W 

Q 1 UQ 2 ,t ,e ^ a Q l ,r,e^ a Q 2 ,r,e • 

Therefore the following limit exists 

h? e = lim -Ug> T£ . 

Moreover, /i^ is non-increasing in e. This proves that the limit Eq.(2.2) exists. 

We now show that it is independent of r. As in the treatment of topological entropy 
for finite dimensional systems, we start by giving an equivalent definition. Given a positive 
number £, we say that two trajectories A 1 and A 2 in A are ^-separated in Q before time T (with 
time-steps r) if there exists an integer j e [0, T/r] for which 

IIMA)-<M^2)IIl~(q) > C- 

We denote by Rq t t (C) the maximum of the cardinalities of the sets of trajectories which are 

pairwise ^-different before time T. Since the coverings in Wq are of diameter less than (, one 
has 

R Q,T,t(() - N Q,T,r,C/2 ■ 

Let now U be a finite covering of A by balls of radius ( in L°°(Q). Then, if two trajectories 
differ on U at some time, their distance is at least (. Therefore 

N Q,T,r,2C ~ R Q,T,r(0 ■ 

These two estimates imply immediately that 

lim h? £ = lim lim -1 lim l °* R Q^ e) . ( 2 . 3 ) 



Lemma 2.3. There are numbers 7 > 0, T > 1 and C < 00 such that for all L, e > 
satisfyingL > Ce~ l and for all A l and A 2 inAforwhich \\A l — A 2 \\- L00 (q ^ < e, wehavefor 
any0<t< C~ l eL - 1, 

ll^(A)-^(^)IU (Q __ 1{t+1)) <re*e, 

withT > 1. 

Proof. See Theorem 5.1, Eq.(5.8) below. In fact, the statement was already shown in [CE2]. 
Using the above estimate one has easily for r' < r, 

Rq T A e ) < R o , Tr( £r " le " 7T ) , 

and 

R QL ,T, T (e) < R QL+ce _ HT+i) ^ Tl {eT-^) . 

It follows now immediately from Eq.(2.3) that lim £ ^ does not depend on r. This concludes 
the proof of Theorem 2.2. 
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3. Upper Bound on the Entropy per Unit Volume 

It does not follow from what was said in the previous section that we have defined a finite 
number. We now give an upper bound. We first observe that from the sub-additivity argument 
we have 

lim — log R n rp (e) = inf — log N n rp 

Therefore, in order to obtain an upper bound, we can fix a T and vary L. 

Let H e denote the e-entropy per unit volume defined in [KT], (see also [CE2] for the 
application to the case at hand). This means that we first define Nq (e) as the minimum number 
of balls of radius e in L°°(Q) needed to cover A\q (the functions on the attractor A restricted 
to Q), and then 

logiV (e) 
H £ = lim Y ■ 

We define the upper dimension per unit volume d up of A by 

d = lim sup 



U ? £ ^ r log(l/£) 

Note that it follows from [CE2] that d is finite. It is an open question to prove that the limit 
exists. However, we have a bound: 

Theorem 3.1. The topological entropy per unit volume is bounded by 
where 7 is the constant appearing in Lemma 2.3. 

Remark. In terms of the variables which scale like the system, we have 7 = M^. 

Proof. Let 1] > be a fixed number. Let e > be small enough such that for all e E (0, e ) 
we have 

For a fixed e e (0, e ) let L v s > be such that for any L > L we have, 



1/ H < 



< V 



We now give an upper bound on Rq l T t (Q . For e G (0, e ) , choose a finite covering U of 

A in L°°(Q L+(T+1)e _, ) by balls of radius eY~ l e~ lT '/2. By [CE2] and the previous discussion, 
we know that we can find such a covering U of cardinality at most 

exp((L + Ce-\T + l)) d ■ ( V + (rf up + r^T + log(2r £ - 1 )))) . 
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Moreover, from the definition of 7, if A l and A 2 belong to the same element of U , it follows 
from Lemma 2.3 that 

SUP \\$Ml) -®t( A 2)k°°(Q L ) < £■ 
0<t<T 

Therefore, 

RQ tTtT (e) < exp((L + C £ - 1 (T+l)) d .(r ? +K p +r ? )( 7 T + log(2r £ - 1 )))). 
It follows that 

lim ~~~t inf -j- log R n ^ ^ (e) < lim -^r -L log R n ^ ^ (e) 
and the result follows by letting T — > 00 and then 77 — > 0. 



4. Localization in Momentum Space and Bounds on the Semi- 
Group 

In this section, we deal with some simple bounds on the semi-group generated by DA. We 
begin by constructing the localization in momentum space. Let x > be a smooth function 
with support in \k\ < 2 and which is equal to 1 for | k | < 1. We also assume x < 1. We shall 
denote Xk<k = an d Xk>k = 1 ~~ x(k/k+). We define the convolution operators 9 r , 

Sr,> and S T) < by 



S T (x) = 


^ dA;e ifcx - 


S Tl >(x) = 


^ d£;e* fcx e 


S T ,<(x) = 


1 dke lkx e 



x . 



Notation. The constants Bq , Sj", . . . , do only depend on the quotient D^jv^ (something like 
the condition number of the matrix D), but not on any other parameters of the problem, except 
d. We also recall that the constants C , . . . , are numerical factors which do not depend on any 
parameters of the problem except d. 

We now state and prove various estimates on these kernels. 



Lemma 4.1. For every p > there is a constant Bq(p) such that for all r > one has the 
bound 



/ 
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Furthermore, convolution with S r is a well-defined operator on bounded functions and has norm 
(as a map from L°° to itself) bounded by 



(4.2) 



for some 



Proof. For simplicity, we write the proof for the case of d = 1, but with a distinction of the 
upper and lower bounds (D^, resp. z/J on the matrix D. For example 



/' 



/ a 7„ „ikx —Dk t 

x dke e 



Thus, in this case, 



/ 



dke lkx e~ Dk T 



= J dk (-id k e lkx )e- Dk2r 
= -i J dke lkx e- Dk2r 2Dkr . 



< \x\~ l J dke~ DklT 2D\k\T 

° { \x\ 2) {vj-y/^vj 



(4.3) 



Using D^/u^ > 1, the generalization to arbitrary d and p is easy and is left to the reader. The 
second assertion follows by taking p = d + 2 and integrating. The reader can also check that 

B\ =C l (Djuy+\ for some C v 

Lemma 4.2. Let r > 0. For all p > there is a constant B% (p) such that for all z G C d one 
has a bound 



I 



dke^e~ Dk r Xk<K 



< Bl{p) 



1 



2k Jim; 



(fc.- 2 + iv-)*/* (i + ^—y /2 



(4.4) 



where \lmz\ = Yli=i \^ mz i\ 
Proof. We get, for z G C, 



zP Jdke ikz e- Dk \ k<K = Jdk{(-zd k re^)e- Dk ^ X (k/K), 
and integrating by parts this leads to a finite sum of terms of the form 



0(1) / dke' kz e- DkT (DTk) ni k-^k-^(d k 3 X )(k/K) 
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where n l +n 2 + n 3 = p, with n 2 < n l . We let / = <9^ 3 x, and we write n l = s { + s 2 , n 2 = s { , 
and n 3 = s 3 , where now s { > and p = 2s 1 + s 2 + s 3 . Thus we need to bound expressions of 
the form 



J\k\<2k t 



(4.5) 



It will be useful to consider first the case v^r < \jk\. Then we can bound (4.5) in d dimensions 
by 

o(i)(£±yki(v.T)'i +a *k?k-'> < o(i)(^)^x 2s '- 2s2 ^r 3 = 0(1)(— ) P kt P - 

^* z/ * ^* 

(4.6) 

In the case when v^r > 1/kl, we bound (4.5) by 

0(1)(— ) p f d/ce-^ fc2r Kr) Sl+S2 |A;| S2 (^T) S3/2 



D 



< 0(l)^_±)P(jy 7 --)-d/2+s 1 +s 2 -S2/2+s 3 /2 



7/. 



(4.7) 



= 0(l)(^)P(^ r )-^/2+P/2. 



Combining (4.6) and (4.7), and observing that \e lhz \ < e 2fc * |Im2:| on the support of Xk<k > we 
conclude the proof of Lemma 4.2. 

Lemma 4.3. For every p > there is a constant (p) such that for all > and aJJ r > 0, 
one has the following bounds: 
i) When v^t > 1/kl one nas 



I 



dke^e~ Dk r Xk>K 



< b;{ p ) 



~v„k t/2 



K-) d / 2 (i + S) p/2 



(4.8) 



ii) When v^t < 1/kl one nas 



dke ik *e~ Dk r Xk>K 



< B;(p)e-"* k * T/1 
1 



+ 



M^{i + ^y' 2 K d {i + ki\xf) 



(4.9) 



P/2 / ' 



Proof. Assume first that v^r <l/kl. In that case, e u * * T < e, and thus it suffices to produce 
a bound without exponential factor. We can write Xk>k = 1 — Xk<k an( ^ § et a bound by 
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combining Eq.(4.1) with Eq.(4.4). This leads to 



dk e lKz e 



kz „ — Dk t „ 



Xk>k t 



< eB*(p) 



1 



+ eB*(p) — 



K^(i + Jgi)^v 
1 l 



D 



(4.10) 



(^ + ^/ 2(l + _l^) 



P W 2 V i/. 



from which the first assertion follows. In the case v^r > l/k 2 , we integrate again by parts and 
get to bound an expression of the form (we work again in the case d = 1 only): 



x I dke^e- Dk ^Xk>k 



= -i jAke^e- Dk ^{2Dk{\- X {k/K))-K l {d kX ){k/K)) . 
Therefore, we get, as in (4.3), 

< \x\- l Q(l)e-"* k2T/2 J dke- Dk2T/2 \2Dkr\ 
+ \x\- l O(l)e-^ 2 [ dke- Dk2r / 2 \K l X '(k/K)\ 



(4.11) 



dke^e~ Dk r Xk>km 



< cr-^ k2T/2 ( D * T V /2 1 ( D *yn 



+ Ce~^ k T/1 



1 k~ l 

(^ryn \ x \ 



< Ce -.* 2 r/2(^iI)l/2 l^/^*)l/2 

The generalization to arbitrary p and d is easy and is left to the reader. 



5. The Dissipative Bound 

In this section, we consider in detail the equation 1 

w(x,t) = DAw(x,t) + M(x,t)w(x,t) , 



(5.1) 



where w takes values in H N . Our bounds will work in dimensions d < 3. We first state the 
assumptions of the Introduction in a more precise form. We first assume that 



||M(M)|| < M,, 



(5.2) 



For simplicity we assume isotropy of the diffusion in the d components of the coordinates, but this requirement 
could be dropped if desired. 
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for all x e R d , t e R + . We assume further that 

\w(x,t)\ < Q # , (5.3) 

for all x, t. (This is the reason for the choice of in the bound on it.) Here, and in the 
remainder of the paper, | • | is the Z 2 -norm of a vector in R^. 

We will fix the constant k^ only in the next section. But we will work here with the 
following "comparisons of scales" which will be essential in the bounds: 

r. = , (5.4) 

Si = , (5.5) 
kl > 1/KrJ . (5.6) 

We consider next the integral representation of w x (•) = «;(•, t = t+): 
(x,tJ = {9 T .*w)(x)+ [ T *ds{9 Tt _ s *{M(.,s)w(;s)))(x) = w[ l) +w[ 2) . (5.7) 

JO 



w 



Using the decomposition 9 T = S r < + S T >, we split w^' and w[ 2 ^ into high and low frequency 
parts: 

(1) , (2) (1) , (2) 

w l,< = W i,< + W l,< > W i,> = w l,> + <> ■ 

Note that w K has Fourier components in {\k\ < 2k^} and w > has components in {\k\ > k^}. 

Theorem 5.1. (Dissipative bound) In dimension d < 3, we have the following bounds for 
some constants B\ and B* 5 : Upon localizing in space, we have for < r < t+, 

sup \w(x,r)\ < B 4 *(sup \w{x,0)\ + VJt)\W;0)\\oo)- (5 8) 

\x\<X~e |x|<A V ' 

Localizing in position space and momentum space, we have 

sup |w > (x,rj| 
\x\<\-i 

< b; ( e -(»A/M*)n + ^y-dA (5 . 9) 
• ( sup \w{xM + {K/m<MU- 

\x\ < X 



Proof. We define, as in [CE1], a family of space cutoff functions: Let 

= l + {{x _ Z a)2/5l y +d/2 = iK* ~ «) • 
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where Z = 5~ d j J dx(l + x 2 )~ 1 ~ d ^ 2 is chosen such that J dxtp(x) = 1. We start with a bound 
in L 2 : 

Lemma 5.2. There is a constant B£ such that the solution of Eq.(5. 1 ) satisfies: 

sup sup / dx^ a (x)\w(x,t)\ 2 < B;( sup \w(x,0)\ 2 +\\w(-,0)fooVM 2 ) ■ (5-10) 

0<t<r„ \a\<\-£ J \x\<\ 

Proof. Let 

X t = j dx^ a (x)\w(x,t)\ 2 . 
Then we have, from the equations of motion, 

d t X t = J dx ip a (w ■ (DAw + Mw) + (DAw + Mw) ■ w) 

= -2 J dx^ a (Vw) ■ D &ym Vw + 2 J dx V> a M syn >| 2 ~ J dx WJVu; • Dw . 

Here, M sym = (M + M T ) /2 and || M || < ||M|| < M^. Observe that by our choice of if) we 
have | < C 2 ip/S^ for some constant C 2 independent of 5^. Using the definitions (1.3) of v % 



and D^, we find 



\d t X t \ < -2u,J dxiP a \Vw\ 2 + 2 j dxif; a Mjw\ 2 + 2C 2 J dx (Dj 8^ a \w -Vw\ . 
We polarize the term containing w ■ Vw and use the identity (D^/5^) 2 /u^ = M^. see Eq.(5.5), 



Then we see that we can find a constant C 3 such that 



\d t X t \ < C 3 M.(l + {DJv.))X t . (5.11) 
Coming back to the assumptions of Lemma 5.2, we see that when \a\ < A — £ we have 

X = / dx^ a \w \ 2 = / dxi(j a \w \ 2 + / dxif; a \w \ 2 

J J\x\<\ J\x\>\ 

< sup | w ( x ,0)| 2 + C 4 ||«; || 2 oo (^) 2 . 

Using this bound on the initial condition, the differential inequality (5.11), and r^M^ = 1, the 
assertion of Lemma 5.2 follows with B£ = exp(C 3 (l + (.D^/z/J)). 

We begin the proof of (5.9). We deal first with the bound on wf'y (x), when \x\ < A — £. 
Consider, for > t > s > 0, the quantity 

Y t s = J dy9 t _ s > {x-y)M(y,s)w(y,s) 

= [ % ^j== V } V iAx - y)M(y, s)w(y, s) . 
J y/ip(x-y) 
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We consider first the case v^t > : Then, by the Schwarz inequality and Eq.(4.8), we have 

Yl < J ^ % ~^ X S y f ^ 1 J dzi; x (z)\w(z,s)\ 2 



< 



f B 

dy- 



■ Ml I dz^ x (z)\w(z,s)\ 2 . 
Since we deal with |x| < A — £, we get from (5.10), 

where K = sup| x | <A \w(x, 0)| and > \\w(-, 0)11^. Note now that by Eq.(5.5), 

Si = > D.(t-s), 

since > t > s > 0. Therefore, 



l+d/2-p 



• M^(K 2 + g 2 ,(^) 2 ). 

Taking p = d + 2, integrating over y, and using again 5 2 = D^/M^, we get, for some constant 
B*\ 

2 



x e -".fc,(*-«) 



(5.12) 



Z/* /2 (t - S ) d / 2 



Taking the square root of this bound, integrating over s, (and using at this point the hypothesis 
d < 3) we get for all r G [0, rj and all |x| < A — 



l«£ 2 W)l < ^ T d.|i; )S | < dsB*Ml- d '\vX) d/ \K + Q^-j) 



e~ u * k * s/2 
t > {vMs) d/A 

(5.13) 
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(2) 9 

This completes the study of the contribution of w\ > to the bound of Eq. (5.9), when u^t > \/k ie . 

We next deal with the case v^r < l/k 2 . The contribution corresponding to the first term 
of Eq.(4.9) is treated as before noting that in Eq.(5. 13) we actually integrate over all r, and that 
the inequality v^t > \/k 2 was not used anywhere except for being able to use (4.8) which is 
the same as the first term in (4.9). Thus, we consider here only the second term. Although the 
bounds are quite similar to the previous case, it might be better to actually spell them out. By 
the Schwarz inequality and Eq.(4.9), we now bound 



,. , B?(F)V-- t -"-"(^) 2 ''l + (|,- i< |V^)'^ 



/' 



• Ml J dziP x (z)\w(z,s)\ 2 . 
Since we still deal with \x\ < A — I, we get from (5.10), 



2 , 



Ks < Bt{p)\?±y»Z-i J dy e ~^ d ~ S) (1 + kl\x- y\r? {l + y +d/2 



Note now that by Eq.(5.5), 

Si = DJM > k~ 2 . 

This leads to 

Ks < (B;(P)) 2 (^ Y P Z~ l [ dykl d (l + kt\x-y\ 2 y +d / 2 -r 

J 

■M 2 B* 6 (K 2 + Ql( S -f) 2 ). 
Taking p = d + 2 and integrating over y we get this time, for some Bg : 

w 2 yS < {B;) 2 5 d K d kl d Ml{K + Qjf) 2 
= {B* 9 ) 2 M 2 (5M d (K + Qj-f) 2 . 



(5.14) 



Taking the square root of this bound, integrating over s, and noting that vJz^t < 1, we get 
r<> s \ W T,,\ < £ dsB;M,(6,k,f 2 (K + Q^) 



B;& d/i (^%y- d,1 (K+Q, s j)- <5 ' l5> 
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This completes the bound for vJs+t < 1. 

The contribution from luj 1 ^ is easier to bound. Using the definition of w^y, we split 

the convolution integral into the region \y\ < A and its complement. Because v^t^ > k~ 2 by 
Eq.(5.6) we can use the bound of Eq.(4.8), and we get for the first contribution when p > d + 1, 



dy 



w(y,0) 



< B* e-"* k * T * /2 sup \w(y,0)\. (5.16) 



For the second term, where \y \ > A, the restriction of the bound to \x\ < A — £ implies \x— y\ > £, 
and then we get using = 1 /M # : 



f 

J\x- 



-v t k T./2 



(Dt) 1/1 ,ne-"* k2 * rJ2 

^ (KT*) d/2 y 1 



(5.17) 



Combining Eqs.(5.13)-(5.17), the inequality (5.9) follows. 

The proof of Eq.(5.8) is very similar to the one given above, and we indicate just the few 
modifications needed. Instead of the kernel % T > we now use the kernel Q T , and we will call Z t 

the quantity corresponding to Y t s but with S T replaced by S r >• Consider first v^r > l/k\. 
Since the bound (4.8) is of the same type as the bound (4.1), but without he exponential factor, 
all bounds go through as before up to the inequality (5.12), which is replaced by a similar one, 
but without the exponential factor. Taking again the square root and integrating over s, we get 
for d < 3, 

< (2 W)| < [ T ds\Z T!S \ < f T *ds(B*) 2 
Jo Jo 



Ml- d/4 (uX) d/4 (K + Q^-f) 



1 



(5.18) 



< B* u Ml- d '*{K + Qjf)Tl- d l A 

= b* 2 (k + qJ-±). 

This bounds the contribution of wf^ to Eq.(5.8) where u^t > In the opposite case, we 

argue exactly as in the proof of (5.9), since the exponential factor in (4.9) was anyway of no 
use before. Finally, the contribution of is bounded in exactly the same way as the one of 



w 



(i) 
L> 



, except for the exponential factor, and we get the bound (5.8). The proof of Theorem 5.1 



is complete. 
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6. The Sampling 

We assume that w l = w(-,t = rj satisfies the "sampling bound" 
where 

|/|^ )A = sup , (6.2) 

RZ" : <5Jn|<A} 

withe < Q^. Furthermore, we assume 

sup \w(x,0)\ < K . (6.3) 

\x\<\ 

Then we have the 

Theorem 6.1. Consider the solutions w(-, t) of Eq.(5.1), and assume that the bounds (5.2), 
(5.3) hold. There are constants B^ 3 , B* 4 , and B^ 5 independent of e, K, and A such that the 
following holds: If the initial condition satisfies (6.3) and iu(-,tJ satisfies Eq.(6.1) with the 
definition Eq.(5.5) of 5^, then 

sup \w(x,t^\ < K' , 

\x\<\' 

where 

K' = B* 3 e + 3K/4 , (6.4) 

and X' = A - B* 4 - B* 5 /K. 

Remark. Since 3/4 < 1, Eq.(6.4) shows that K' is smaller than K as long as K > AB^e. 
(In fact, by choosing different constants — in particular a large — we can achieve any ratio 
p > instead of the 3/4.) 

A corollary of the proof is the following result which can be viewed as a generalization of 
the dissipative bounds of Theorem 5.1. 

Corollary 6.2. In dimension d < 3 we have the following bounds 



3 

sup \w(x,rj\ < B* 16 \w(-,tJ\ s x + - sup \w(x,0)\ 
\ x \<\-e 4 |x|<a 

+ B 1 * 7 ^||«,(.,0)|| oo . 



(6.5) 



Proof. Our proof is based on sampling theory for functions in the Bernstein classes, (see, e.g., 
[B]). 

Definition. We call the Bernstein class of entire analytic functions / bounded by 

\f(z)\ < Se ff|Imz| . 
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Our first observation is that S T)< maps into a Bernstein class H> B lk : 
Lemma 6.3. The convolution operator with kernel 

%,<{*- V) = J dke ik ^e~ Dk \ k<K 

is bounded from L°° to itself. Furthermore, for r < r^, it is bounded from L°° to the Bernstein 
class 23 B , 2k ' an d S T < * / is an analytic function of z and one has the bound in C d : 

|(S T>< */)(*)| < B^^WfW^ . (6.6) 

Proof. This follows at once by taking p=d+2 in Lemma 4.2 and integrating. 

Using this result, we can now bound the low-frequency part w 1< of w 1 (x) = w(x, rj: 

Lemma 6.4. There is a constant B* g such that the term w l < is bounded by 

IK<(*)lloo < B* 20 e 2k *^\\w(.,0)\\oc < B* l9 e 2k *^ . (6.7) 
Proof. We recall the representation 

W l,<(-) = W<(;T*) = 9 T „,< *«><) + / * S Tm - a , < *(M(; S )w(;s)) . 

Jo 

By Lemma 6.3, this is bounded as follows: 

!«;<(*, r.)| < e 2k ^{B* 8 Q,+BZ(p) r.M.Q.) = e 2k *^B* l9 , 

so that Eq.(6.7) is proved. 

Now that we have established that the function w l < is entire analytic and exponentially 
bounded, we can use the following sampling result (written for functions in d = 1): 

Theorem 6.5. [B] Assume f £ 23 Scr . The following representations hold: 
/-(0) 4a f, (-1)" 



/(x) = /'(0) sin(^) + /(0)^^ + ^sin(^) £ w Jjl^ /fen) • (6-9) 



ffere r = — 

Remark. In higher dimensions, the sum is over a lattice, and details are left to the reader. 
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Proof. [B]: Theorem 11.5.10 and Eq.(11.3.1). (There is an obvious dimensional misprint in 
11.5.10.) 

Note that the sums in Eqs.(6.8) and (6.9) are absolutely convergent. We now bound them 
as follows. Let N be a (large) integer* . Then, in the first expression, the sum over the terms with 
\n\ > N is bounded by 0(N~ )\\f\\ oc> - In the second expression, we assume for the moment 
that x E (x , x 2 ). Then the sum over \n\ > N is bounded again by 0(N~ l ) 



Proposition 6.6. Assume that f e S S (J . Assume furthermore that f satisfies the bounds 

\f{x j )\<a for \j\ < J, 
where Xj = irj /(2a). Assume furthermore that WfW^ < S' . Then, one has the bound 

\f{x)\ < C 6 a+ for all \x\ < 7r( J — N — 1)/ (2a) . 

The constants C 6 and C 7 are independent of the parameters of the problem. We shall assume 
C 6 >\. 

Proof. We first bound f'(0) using Eq.(6.8), and this leads to 

|n|<JV V ' \n\>N V ' 

If \x\ < Jn/(2a), then x E (x lfl ,x lilJr7 ) where h is the integer part of ax/n. Note that by 
assumption we know that \f(x 2n )\ < ot for all n satisfying \n — n\ < N. Thus, splitting the 
sum as before and shifting the origin to fin /a, we get a bound 

\f(x-(nn/a))\ < V - 1 + ^S' V - 1 

1 v v ' ' n ~ 7r 2 ^ (2n+l) 2 7r 2 ^ (2n+l) 2 

|n|<JV V ' \n\>N V ' 

Eax sin(ax) v-^ 1 

— 7 " — ~~\ + S / • 

\n\<N v ' \n\>N 

The assertion follows. 

We can complete now the proof of Theorem 6.1. We first bound w l > . Setting X" = X — £ 
and using the bound Eq.(5.9), we get 

sup \w h> (x)\ < B; (e-Ml/".)/2 + (^yM 

\x\<\" \ V*K J ^ 

• ( sup |«;(x,0)| + (^)||«;(.,0)|| oo ). 

\x\ < A 



In this section, N denotes just an integer, and not the number of components of u. 
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We now begin fixing the constants: Recall that Q^, and <5 + = D^t^ are given by the 
parameters of the problem. We have also assumed, as a hypothesis of Theorem 6.1, that 

sup |tu(a;,0)| < K . 

\x\<\ 

We choose £ = S^Q^/K, so that the last factor in (6.10) is bounded by 2K. We next choose 
so large that 

B; ( e -(»X/M*)/2 + (^_y- d /4\ < 1 (6 n) 

More precisely, we let 

K = C % \6C 6 (BtY^- d ' 2 \Mj V ^/ 2 = B* 2l (MjvS /2 , (6-12) 

and this is our final choice for k^ . Clearly, if C 8 is sufficiently large, both terms in the sum 
(6.11) will contribute less than 1/(16C 6 ). 

Remark. It is at this point crucial that the construction of the quantities Bq , B\ , . . . , did not 
depend on k^, since all these constants depend — as we have said before — only on the quotient 
D^/v^. In particular, for problems where u had only one component, this would mean that the 
B* are just pure numerical factors, since then =v if . It is also important to note that k^ and 
5^ do not depend on the quantities A or K which occur in Theorem 6.1. 
Thus, so far, with our choices we conclude from Eq.(6.10) that 

sup \w h> (x)\ < -±-K. (6.13) 

\x\<\" 

We next observe that w x < = w x — w x > and so (6.13) leads to 

K<«)| < \wMK)\ + K>(™UI < \w x {n5.)\ + -^-K , (6.14) 
provided 5Jn\ < X" . 

By construction, w l < has a Fourier transform with support in \ k\ < 2k^ and furthermore, 
by Lemma 6.4, it is in S s a with S = B\ 9 and a = 2k^. 

We can now apply Proposition 6.6 to the function / = w 1 <5 with S' = S = B\ 9 . 
Choosing J = \2a\" /ir] (here, [•] is the integer part) we conclude that, for \x\ < X' = 
A" - 6+ - tt(N + l)/(2fc # ), one has abound 

\w h< (x)\ < C 6 \w lt< \ SmtX „ + C 7 N- l Bt 9 . (6.15) 

It is useful to introduce B 22 = C 7 B^ 9 Using Eq.(6.14), this leads to 

sup \w h< (x)\ < C 6 (\ Wl \ SmtX „ + -^K) + N-'B* 2 . (6.16) 

\x\<\' 4U 6 
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We next choose 

N=Ub* 22 . (6.17) 

Then we get 

sup \w lt< (x)\ < C 6 \ Wl \ s x „+K/2. (6.18) 

\x\<\' 

Combining Eqs.(6.18) and (6.13), we get 

K K 

sup \w x {x)\ < sup \w l < (x)\+ sup \w l (x)\ < C 6 \ Wl \ s y + — + — . 

\x\<\' \x\<\' \x\<\" L 4C/ 6 

Since we have taken C 6 > 1, we see that we get finally 

sup Iw^x)] < C 6 \ Wl \ 5 y + 3K/4 , 

\x\<\' 

and 

A' = A" -6* -tt(N + 1)/(2K) 

= A - 5*B* 2 /K -5,- n(N + l)/(2fcj . 

From the definition of N = 0(1 /K) it follows at once that there are constants B\ A and B\ 5 
such that 

V > A - B* u - B* l5 /K . 
The proof of Theorem 6.1 is complete. 

Theorem 6.7. (Sampling bound) Consider the solutions w(-, t) ofEq.(5.1), and assume that 
the bounds (5.2), (5.3) hold. There are constants B^, B 24 and B 25 such that if 

Mi6„t-jrJ\ < £, (6.19) 

for all \i\ <(L + B^/e)/5^ and all \j\ < B^\og{\/e), then 

sup \w(x,t)\ < B 25 e . (6.20) 

\x\<L 

Remark. It will be seen from the proof that it suffices to sample on a somewhat smaller, 
non-rectangular domain for (6.20) to hold. See also Fig. 1. 

Proof. We obtain the result of Theorem 6.7 by iteration of Theorem 6.1. With e given as in 
the statement of Theorem 6.7, let m be the smallest integer for which 

(3/4)"U < e , 

and note that then 

m < log(Q,/£)/log(4/3) + 1 < C 9 log(l/e) . 
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when e < \. The number m will be the number of time steps needed to achieve a precision 
B^e in Eq.(6.20), where we define = 2£? 1 * 3 , and B^ A = C 9 . Let now K = and 
K j+l = B* 3 e + and let L m = L, and L j _ l = L- + <y„ + B\^jK-_ x . Assume now 

H*<M-jrJ| < e, (6.21) 

for all | z | < L m _ i /5 >¥ and for all j = 0, . . . , m. The reader can check easily that our definitions 
are made such that Theorem 6.1 applies at each time step considered. Since we have 

IM.f-jOlloo < Q*> 

for all j, we can inductively bound, for j = 1, . . . , m, 

sup \w{x,t- (m- j)rj\ < C 6 e + 3K J _ 1 /4 = JT. +1 . 

Note now that K- = Q^p J + C 6 e(l - p J )/(l - p), where p = 3/4 and thus 

i=l 

Clearly, there is a C 10 > such that if > C w pP , so that we get 

°10 i=l ^ 

In particular, we can find some Cj t , so that 



L < L + mS, + -^-^—- < L + C 9 dJog(l/e) + C n /e. 
°io P ~~ 1 

Note that L is the width of the "earliest" bound in (6.19), and we choose B^ = C 9 5^ + C n . 
Since K m < B^e we have shown that the bounds of Eq.(6.21) are sufficient to ensure 

sup \w(x,t)\ < 2B* 3 e . (6.22) 

\x\<L 



The proof of Theorem 6.7 is completed. 
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7. Outlook 

In this section, we wish to discuss potential experimental aspects of our results. These aspects 
must necessarily rely on a number of conjectures about the system under consideration, and are 
similar in spirit to the discussion found in [ER]. While we have defined topological entropy in 
Section 2, we now need to address the question of entropy relative to an invariant measure \i. 
We shall call it h . More precisely, let / be a continuous map of a compact metric space and let 
ji be an /-invariant non-atomic ergodic measure. The entropy h (f) is determined as follows. 
For e > and an integer n > let 

V(x,e,n) = {y\d(f i (y),f i (x))<e, < % < n) . 

One has (see [BK], [Y]): 

h (f) = — limliminf — \ogfi(V(x,e,n)) . 

This is not a very convenient expression for explicit computations. As for the case of finite 
dimension (see [GP]) one can try to determine instead a correlation entropy. In our case we 
would also like to include the dependence on the size of the window in which the system is 
observed. One is naturally lead to the following definition of correlation entropy per unit length 
K 2 (n) — and unit time — (see [GP, ER]) for a measure fi which is space and time invariant and 
ergodic. The definition of this quantity is: 

K-y(ii) = lim lim — 7 lim — 

/ 1 N-l n-l , N \ (7.1) 

l0§ U^oo W £ II e ( 6 - ll*W^) - ^(i + k)r(A)\\^( Ql 
\ j,k=0 i=0 V 

where © is the Heaviside function, and <£> t is the flow of the dynamics. 

It is an open question to prove that the limits in the above expression exist for \x almost 
every A (except the limit over e since the quantity is increasing). Furthermore, the result should 
also be independent of the time step r. But let us assume that these limits exist. In that case our 
sampling bound Theorem 6.7 gives us a constructive handle on computing the r.h.s of Eq.(7.1). 
We recall the definition 

|/|, >A = sup \f(n8)\. (7.2) 

{nEZ d : S\n\<\} 




Theorem 7.1. Assuming the limits in (7.1) exist for ^-almost every A, we have for every 
5 e (0, 5 J and every t G (0, rj: 

K (n) = — — lim lim -^j lim — 

/ AT-ln-1 , ,\ (7.3) 

\ j,k=0 t=0 V 7 / 
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Remark. Note that the sum 

Efjio nr=o ©(•••) counts the number of pairs of points in the 
sample which have a "distance" of less than e in embedding dimension n, where the distance 
is measured with the discrete sampling step 5 over a region \x\ < L. In this sense, our result 
says that for PDE's on the infinite line, measuring the K 2 entropy can be done by the usual 
Grassberger-Procaccia algorithm [GP, ER]. 

Remark. If we choose 5 = 5^ and r = r^, then the "sampling error" is controlled by the 
constants of the proof of Theorem 6.7. But if we choose, for example r = rjl, the number 
of time steps — i.e., the factor B 2it — needed in Theorem 6.7 will be the same as if had been 
replaced by 2M^ in the original assumptions. So clearly, the total time during which one must 
measure to achieve a given sampling precision can not be shortened by sampling at shorter 
intervals. Similar precautions are necessary for sampling in space. 

Proof. In view of the remark, we give the proof only for the case r = r + , 5 = 6^, in which 
case the constants B^, . . . , retain their meaning from earlier parts of the paper. We obviously 
have 

N-1 n-1 



^ iV — 1 Tl— 1 / \ 

j,fc=0 i=0 V ' 
. N-1 n-1 , x 

>^E II i £ -\\^)rM)-^i + k)rM)\\^(Q L )) ■ 
»■ i— n i— n \ ' 



j,k=0 i=0 

On the other hand, applying Theorem 6.7, we have for n > B 24 log(l/e) + 1 and L > B^/e: 

n-1 



i=0 

n-1 



n e ( £ -i^),^-w,^)ii,L) 

i— n \ / 



< n e ( s 2 v - ii<w.w - *( l+fc K(^)ii^(Q L _ fl , /£ ) 



\ 23 

^ II e f S ^ - ll*(i+i+[Bllog(l/ e )])r.(^) 



=[B 2 * 4 log(l/e)l 
n-l-[B 2 * 4 log(l/e)] 



i=0 



(A) ||, 



(i+k+lB^ log(l/e)l)r. W \\L.°°(Q L - B £/e) 



Therefore, 

N-1 n-1 



^ i\ — i n — i s \ 

^2 E II e e - |$ ^ (A) "V)r.^lw 
t—n i— n \ / 



j,fc=0 i=0 

7V_1 n-l-[B 2 ;iog(l/e)] 



%E II e(Bie-||* (4+i)T .(^)-* (4+fc)T ,(A)|| L „ 



j,fc=0 i=0 

+ lo(l+log(l/ £ )) 2 . 



'(<Jt-fl'/«) 
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The result now follows by taking the limits N — > oo, n — > oo, L — > oo (in that order) for e 
fixed, and then letting e tend to zero. 
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